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We derive a microscopic model for dissipative dynamics in a system of mutually interacting 
qubits coupled to a thermal bath that generalises the dissipative model of Landau-Lifschitz-Gilbert 
to the case of anisotropic bath couplings. We show that the dissipation acts to bias the quantum 
trajectories towards a reduced phase space. This model applies to a system of superconducting flux 
qubits whose coupling to the environment is necessarily anisotropic. We study the model in the 
context of the D-Wave computing device and show that the form of environmental coupling in this 
case produces dynamics that are closely related to several models proposed on phenomenological 
grounds. 


I. INTRODUCTION 

The Landau Lifschitz-Gilbert-equation^’^ provides a 
phenomenologically motivated description of the stochas¬ 
tic, dissipative evolution of a spin system. Conceived as 
a model for an open magnetic system, the dynamics con¬ 
sists of two terms corresponding to processing, Hamilto¬ 
nian evolution, and noisy, relaxing, dissipative evolution. 
This form of dissipation corresponds to the classical limit 
of a spin system coupled isotropically to a bosonic bath. 

The behaviour of a Superconducting quantum device 
can be mapped to spin dynamics^. In the case of a flux 
qubit, the low energy dynamics map on to those of a 
spin-1/2 magnetic moment, and an isolated qubit will 
exhibit only the Larmor precession described by a dissi¬ 
pation free Landau-Lifschitz-Gilbert equation; the qubit 
wavefunction evolves such that the absolute amplitudes 
of the two energy eigenstates remain constant whilst the 
phase between them changes at a constant rate. 

This relationship suggests that qubits coupled to the 
environment may display the same dissipative behaviour 
as magnetic moments, governed by the full dissipative 
Landau-Lifschitz-Gilbert equation. This model has been 
used to model an extended array of superconducting 
qubits in Ref 4 whilst related vector models have been 
used in Refs 5-7. 

However, the Landau-Lifschitz-Gilbert equation de¬ 
scribes the dissipative dynamics of a two-level system 
that is exposed to isotropic environmental coupling, i.e. 
identical baths coupled to the s^, Sy, and operators. 
For qubits, these operators may have different physical 
origins and hence will couple differently to noise. Thus, 
the dissipation will also be anisotropic, as the stochastic 
noise and dissipative terms are related by the fluctua¬ 
tion dissipation relation. As a result, systems, such as 
flux qubits, with anisotropic couplings have correspond¬ 
ing anisotropies in the dissipation and noise. 

Due to the physical geometry of the superconduct¬ 
ing flux qubit, stray flux and other environmental effects 
couple to the Sz operator. The anisotropy of this envi¬ 
ronmental coupling introduces qualitatively new features 
into the system’s dynamics. We show the existence of a 
regime where the qubit dynamics are typically confined to 


a low dimensional sub-manifold of the full Hilbert space. 

The dissipative dynamics of a two level system has 
been studied extensively^. In this work, we use a 
Keldysh path integral®^^® to derive a Langevin^^’^^ de¬ 
scription of the dynamics of a flux qubit accounting for 
the anisotropic coupling to the environment. This ex¬ 
tends the previous work on dissipative spin models^^^^^. 

An accurate model for the dissipative dynamics of a 
flux qubit can be used to assess the capabilities of puta¬ 
tive quantum technologies. We apply our model to the D- 
Wave computing machine, which consists of a large array 
of controllable flux qubits. Extensive analysis has sought 
to correlate the behaviour of this machine with various 
quantum and classical models^“^’^^“^®. Since classical dy¬ 
namics correspond to a particular restriction upon fully 
quantum dynamics, the effectiveness of this approach is 
dependent upon identifying the appropriate restrictions 
that correspond to the classical limit. We show in an 
appropriate strong coupling limit the biasing of trajecto¬ 
ries in the anisotropic Langevin equation allows one to 
obtain dynamics reminiscent of the heuristic models of 
Refs 5, 7, and 24. 


II. THE ANISOTROPIC DISSIPATION OF 
FLUX QUBITS 

A superconducting flux qubit will couple to its envi¬ 
ronment in various ways. The environmental degrees of 
freedom may consist of charge fluctuations, flux noise, 
coupling to nearby spins, or trapped vortices. When it 
is not possible to interrogate such environmental sub¬ 
systems, the information contained in the state of these 
subsystems is lost from the system of interest—the qubit. 
This decoherence inhibits the ability of a qubit to remain 
in a given state indefinitely with good fidelity, and in turn 
limits the ability of an array of qubits to sustain entan¬ 
glement. 

This loss of unitary evolution generates dissipative 
dynamics that, over sufficiently long times, drive the 
system towards certain equilibrium states, or dynami¬ 
cal fixed points. In this sense, the effects of decoherence 
are inherently inhomogeneous over the system’s Hilbert 
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space, they affect different states in different ways, driv¬ 
ing them towards different fixed points or along different 
paths to such fixed points. 

We anticipate that a system exchanging energy with 
an environmental bath will relax to a fixed point given 
by the Gibbs state. There remain, however, a plurality of 
dynamics that result in the system relaxing to this state, 
and we expect that different baths and couplings with dif¬ 
ferent physical origins will result in different relaxation 
dynamics. This is relevant to controlled quantum sys¬ 
tems in which the dynamics of interest occur before the 
system has relaxed to thermal equilibrium, but not nec¬ 
essarily on timescales where dissipation can be neglected. 

The Landau-Lifschitz-Gilbert describes the dissipa¬ 
tive dynamics of a two-level system with isotropic cou¬ 
pling to the environment. It is often quoted in one of two 
equivalent forms; 

s -I- s X [(B -I- 77 ) - 7s] = 0, (la) 
^ [(B+ » 7 )+ 7 S X (B-hr?)] = 0. (lb) 

These are non-linear differential equations in a vector s 
of magnitude s that parametrises a qubit spin coherent 
state s|s) = s|s). The equation is stated here in terms of 
an effective magnetic field, B. More generally, the mag¬ 
netic field may be replaced by an appropriate derivative 
of the spin Hamiltonian, B = — Vsi?(s). rj describes a 
stochastic noise that satisfies the fluctuation dissipation 
relation; {r]a{(jj)r]i 3 {uj')) = yw coth (^) 5ap5{u) + w')^”^. 

This description closely resembles the phenomenologi¬ 
cal Bloch equations used to model nuclear magnetization, 
the difference being the Landau-Lifschitz-Gilbert equa¬ 
tion describes a single spin and not an ensemble average, 
thus the decay must preserve the spin |s| = s. Neglecting 
zero-point fluctuations, this Langevin equation describes 
the dissipative evolution of a single qubit in the presence 
of a magnetic field, B, and coupled isotropically to its 
environment. 

The dynamics of a system of multiple interacting 
qubits may be described by a set of coupled Landau- 
Lifschitz-Gilbert equations. These dynamics are general 
and constitute a restriction to product states. Coher¬ 
ent dynamics are permitted for individual spins in these 
states, but there is no entanglement between spins. In 
this sense the equations correspond to the classical limit 
of the system. 

Although originally developed to describe spins that 
are isotropically susceptible to noise, the Landau- 
Lifschitz-Gilbert equations have been used to model the 
classical dynamics of dissipative qubits^, whilst other au¬ 
thors have also made use of similar vector models^’^^. 
However, as argued above, the physics underlying super¬ 
conducting flux qubits implies that noise and dissipation 
will be anisotropic. Whilst the microscopic origins of flux 
qubit decoherence are not fully understood^®, many en¬ 
vironmental interactions are often modelled using linear 
Sz couplings.® These include flux noise and coupling to 
impurity spin^®^®® or boson®^’®® degrees of freedom, and 


have been shown to be significant contributions to deco- 
herence®"^. Whilst in the context of the D-Wave comput¬ 
ing machine it has been argued that linear coupling is 
the correct minimal model®®. Without loss of generality 
any linear coupling can be chosen to be Sz and whether 
longitudinal or transverse. Due to the close relationship 
between the damping and noise terms enforced by the 
fluctuation dissipation relation, this cannot be remedied 
solely by an appropriate adjustment to the noise term; 
an adjustment to the noise term effects the damping to 
term in a manner that is difficult to guess. 


III. A DYNAMICAL EQUATION WITH 
ANISOTROPIC DISSIPATION 


Here we follow the approach of Ref 8 to find the par¬ 
tition function Z = f for an open quan¬ 

tum system, and evaluate the dynamics of a qubit using 
a stationary phase approximation. In a closed system, 
the action of the spins is given by 


Sn = / dtLn = / dt 




Si 




( 2 ) 


where A^ = ^ is the single monopole vector 

potential. This action describes the precession of a 
spin around the axis of the magnetic field. When the 
spin is coupled to a bath of oscillators through the s^- 
component, the action becomes 


S = So+ J2 / df 


/ • 2 2 2 \ 


(3) 


We use Keldysh field theory to find the dynamics of this 
system. This provides a methodology for treating an 
open quantum system, enforcing the necessary fluctua¬ 
tion dissipation relation between the stochastic noise and 
deterministic damping terms induced by decoherence. 

A state vector is sufficient to define the behaviour of 
a closed system. Non-equilibrium, open quantum sys¬ 
tems require an ensemble of state vectors, or a density 
matrix, because of the loss of information to the bath. 
Evolving a density matrix p requires both pre- and post¬ 
multiplication by the time evolution operator U{ti,t 2 ) = 


T exp 


iSlldtH{t) 


. Thus pt 


U{t,0)poU{0,t) in con¬ 


trast to the time evolution of a state vector which requires 
only one time evolution operator jf/jt) = U{t, 0 )|'i/'o)- This 
gives rise to a doubling of the degrees of freedom in the 
Keldysh theory of open systems compared to those re¬ 
quired to describe the evolution of a closed system. These 
degrees of freedom correspond to the forward and back¬ 
ward branches of the Keldysh contour, C, over which 
Eqs (2) and (3) are integrated. 

It is useful in calculations to separate diagonal and 
off-diagonal contributions to the density matrix. The for¬ 
mer are described by the sum of fields on the forwards and 
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backwards Keldysh contour, and the latter by their dif¬ 
ference. These helds are usually called the classical and 
quantum fields respectively. Expanding in the quantum 
fields generates a series of higher order quantum correc¬ 


tions to the classical dynamics. Performing this rotation, 
and expanding to first order in the quantum fields results 
in each spin having a total action S = Sq + S'diss 



d dLp 
dt dsi 




[D-YJ 



■ (4) 


The first term in Eq (4) encodes the Hamiltonian (closed 
system) dynamics of each spin. The remaining terms 
encode the dissipative dynamics induced by interactions 
with the bath. 


A. A Langevin Equation 


The dynamical equation can be obtained from a limit 
of the Keldysh Theory as follows: In an open system the 
state of bath cannot be interrogated, thus it is integrated 
out to obtain the dissipative contribution to the action 
in terms of the spin alone; 


*5'diss 





0 

DR ijK 



(5) 

where the correlators Qi similarly for 

and , are equal for each spin, and defined by the 
bath spectral function J(w) = a/'^ ^ ~ 

The equation is made linear in the spin fluctuation vari¬ 
ables s* by a Hubbard-Stratonovich transformation 


Sdiss = ^ dt(2 - ?7j [D^] ^ 

- jm. 

( 6 ) 

which introduces the variable ry that later appears as the 
stochastic term in the dynamical equations, s'^ now act 
as Lagrange multipliers and can be integrated out of the 
partition function to obtain 


6 X s^ + -I- J dVD^{t - t')si^z{V)z + rjizJ 


where the Keldysh correlator {t — t') = 

—2i{rii{t)r]i{t')) defines the Gaussian noise correla¬ 
tions. Evaluating this via the fluctuation dissipation 
relation = {D^ — D^) coth (^) one arrives at 

(77,(t)r7j (t')) = ^ J dwJ(u;)coth(^) cos(a;(t-t')), 

(9) 

whilst the retarded correlator defines the damping kernel 

/ t -| pOO 1 

dt'D^{t') = — / — J(w) cos(wt). (10) 

-oo Jo ^ 

In Eq (8) the bare capacitance has been renormalised by 
a term generated by the bath in the usual way^. 

Many of the dynamical features of this model can be 
seen within the Markovian approximation in which the 
bath memory is neglected. In this regime the dynamics 
are described by the Langevin equation which can be 
expressed in two equivalent forms, where = ry^z and 
Bi = —Vs^iL = hi -I- ^ijSj -|- ..., we obtain 

Si -t Si X [(Bj -I- r]i) - 7 z(z • Si)] = 0, (11a) 
Si -I- Sj X (Bi -I- ? 7 j) -1- 7 Si X z [z • (si x Bi)] = 0. (lib) 

The Langevin equation, Eq (11a), is a stochastic differ¬ 
ential equation whose solution is given by an ensemble of 
pure state trajectories. While it is possible to obtain the 
many body dynamics within any manifold of quantum 
states, to obtain the classical dynamics of a many qubit 
system, we have restricted the pure state trajectories to 
product states, thus encoding only classical correlations. 

The Langevin equations, Eqs (11a) and (lib), ob¬ 
tained by this method include, by construction, the dis¬ 
sipative and stochastic effects induced by an anisotropic 
bath. These dynamics are different from those obtained 
for an isotropic bath,^® Eqs (la) and (lb), with impor¬ 
tant consequences. 


the dynamical equations describing the semi-classical dy¬ 
namics of the system are then found are then found by IM. DYNAMICS OF THE MODEL 

taking the saddle points of the path integral. The result¬ 
ing dynamical equation is We now explore the novel dynamics in the presence 

of anisotropic coupling to the bath in its different pa¬ 
rameter regimes. This model, Eqs (Ha) and (Hb), de¬ 
scribes the dissipative dynamics of a system of interact- 


s, = s, X 


V, JL 


dt'r(t — + rjz 


( 8 ) 
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ing, non-entangled flux qubits, with environmental cou¬ 
pling solely through the Sz operator. There are important 
differences between the effects of isotropic environmental 
couplings, Eqs (la) and (lb), and anisotropic couplings, 
Eqs (11a) and (11b). The energy conserving dynamics of 
these two models are the same, consisting of precession 
about the external held, and they both relax to the same 
thermal equilibrium distribution. We shall concentrate 
upon the case where thermal fluctuations are small in the 
sense that the equilibrium thermal distribution subtends 
a small solid angle on the Bloch sphere. This requires 
that (d^) PiT/B TT^, or alternatively B ^ T. In this 
limit, we may sensibly discuss a large deviation from ther¬ 
mal equilibrium and consider the dissipative relaxation to 
it. These dynamics are very different in the pressence of 
an isotropic, Eqs (la) and (lb), and anisotropic coupling, 
Eqs (11a) and (11b), to the bath. 

The Landau-Lifschitz-Gilbert equation constitute a 
special case where the Hamiltonian and dissipative dy¬ 
namics separate in the equations of motion. This sep¬ 
aration is clear when written in polar coordinates. Eor 
brevity we neglect stochastic terms and let B || z, we find 


'jsB 
1 -I- 


sind. 


B 


1 -I- 7^5 


2 „2 ■ 


( 12 ) 


These equations generate motion in perpendicular direc¬ 
tions and, there is no interplay between their dynamics. 

In contrast, in more general dissipative models there 
is no such separation of the effects of precession and dis¬ 
sipation, and their interplay remains important. In the 
Langevin equation for Sz coupling, eqn. (11b), the sys¬ 
tem relaxes indirectly, through the interplay of dynamics 
and the state dependent modulation of the rate of dissi¬ 
pation. The effect of this interplay is highlighted by the 
appearance of regions of novel behaviour, visible in fig 1 , 
such as dissipation free precession, retrograde motion, 
and effective dimensional reduction. 

To discuss the dynamics of the Sz-coupling model, it 
is useful to introduce the timescales = B — |B| and 
= 'jsB sin^ 9* where 9* is the polar angle of the field 
B from z. Tp and Td which are characteristic of the pro¬ 
cessional motion and dissipative motion, respectively. In 
the limits where these scales are widely separated, the 
system’s behaviour is dominated by the faster dynam¬ 
ics on short timescales, whilst some effective dynamics 
emerge on longer timescales. 


A. Weak coupling limit, Tp <C Td 

For weak coupling, the system’s behaviour remains 
dominated by the precession found in the fully closed 
system dynamics. As shown in fig. la, the spin will gen¬ 
eral perform many rotations about the magnetic field, B, 
before reaching the neighbourhood of the ground state. 
This allows the dissipation and noise to be averaged 
over these rotations. In this regime the Langevin equa¬ 
tion may be approximated by a Landau-Lifschitz-Gilbert 


z z 




FIG. 1. Dynamics of spins in the presence of anisotropic dis¬ 
sipation. The direction and rate of evolution over the sphere 
surface are indicated by the streamlines and grayscale den¬ 
sity, darker regions indicate faster evolution, (a) For weak 
anisotropic coupling the state precesses similarly to isotropic 
coupling, (b) For strong anisotropic coupling dynamics is 
markedly different. The system rapidly relaxes to a reduced 
0(2) manifold where it undergoes constrained dynamics. 


equation 

s -h S X [(B -I- ?7eff) -f 7eff s X (B -f r7eff)] = 0 (13) 

where 7 eff = 57 sin 0* is the effective dissipation, This 
results in dynamics characteristically similar to the 
isotropic case (Eqs (la) and (lb)), thus the limit Tp <C Td 
offers no novel dynamics. 


B. Strong coupling limit, Tp ';$> Td 

In the limit Tp ^ Td, when the effects of anisotropy 
show up most profoundly, the system’s dynamics are 
dominated by the dissipative term. This term is state 
dependent and drives rotation about the z axis, however 
the dissipation goes to zero when B, s and z are coplanar. 

Exactly how this reveals itself in the dynamics de¬ 
pends upon the bath memory and temperature, though 
the net effect is similar in all cases: we see a separation 
of timescales and a reduction of the full 0(3) qubit dy¬ 
namics to effective 0 ( 2 ) dynamics characterised by fast 
decay towards, or oscillation around, a reduced manifold. 

The Markovian approximation: In the first in¬ 
stance, it is easiest to analyse our model within an Marko¬ 
vian approximation. Separating Eq (11b) into the slow 
9 and fast </> dynamics we obtain 

9 = —Bsm9* sincf), 

(j) = B (cos 9* — sin 9* (cos (p cot 9 sj sin 9 sin (p)) 4 - 77 . 

(14) 

where p typically relaxes with a characteristic timescale 

~ 7sH Sind Sind* 
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during which time 9 makes a very small change. Thus 
we can treat the appearance of </> in the 9 dynam¬ 
ics as a stochastic variable sampling the quasi-static 
distribution of the (j) dynamics. This distribution is 
given approximately by p{4>) oc exp [—^(0) cos(^ — (/)*)] 
where A{9)~^ = yTr^ cost/)* and tan(/)* = sin(0 — 
0*)/(sysin^ 0sin6**). Equivalent to quasi thermal- 
equilibrium fluctuations in <j). General effective dynamics 
in the slow 9 variable are obtained in appendix A, the 
high and low T cases are discussed here: 

For T B the system quickly relaxes to a state in 
which the (j) distribution is sharply peaked around the 
dynamical fixed point close to (/)* « 0. This confining 
behaviour, as shown in fig. lb occurs when the dissipative 
dynamics drive the system towards a one dimensional 
manifold on the fast timescale Td, after which a much 
slower interplay between the Hamiltonian and dissipative 
dynamics sees the system relax to its ground state. These 
latter dynamics are described by 


Bsin(9*-9) , 

sy sm 9 


(16) 


where the noise {ri{t)) = 0 has correlations {rj'{t)rj'{t')) = 
2jT{T(j,B sin 9*)^S{t — t'). In this limit the stochastic ef¬ 
fect of the bath is weak (in the sense that thermal fluctu¬ 
ations subtend only a small angle on the Bloch sphere), 
whereas it strongly biases the trajectories to dissipate 
energy. In this manner the effect is akin to the trajec¬ 
tory ensemble approach^^ with a transition to dynamics 
confined to </> « ^*. Due to the known relationship be¬ 
tween these dyanmics^®^'^® this Langevin equation can 
also be related to a Monte Carlo 0(2) model on appro¬ 
priate timescales. 

For T ^ B the noise is sufficiently strong that (p 
makes large excursions away from the dynamical fixed 
point. However, due to the separation of timescales, in 
both cases the long time dynamics are captured by an 
effective theory in the slow variable 9. This constitutes a 
dissipative reduction of the phase space from 0(3) to ef¬ 
fective 0(2) dynamics. In the T ^ B limit, this reduced 
dynamics is described by 


sin(6»* - 9) sin6»* 
2Ty sine 


(17) 


where the noise has correlations iv'itWit')) = 
t^{B sin 9*)“^6{t — t'). 

The limit in which the dynamics become Markovian 
is subtle. Often Markovian dynamics can be obtained 
by assuming an Ohmic bath J(w) = 4yw. When the 
temperature is much greater than the characteristic fre¬ 
quency of the system (the frequency at which the bath 
dynamics are sampled) the bath falls in its classical limit, 
coth(w/2T) —>• 2T/w, so that Eq.(9) yields a (5-correlated 
noise. In this regime, however, thermal motion gener¬ 
ally dominates-in the sense that thermal fluctuations 
cover the entire Bloch sphere; yj{9“^)' ~ tt. There are 
alternative-more physically realistic-bath distributions 
for which the Markovian limit arises naturally. 


z 


z 



FIG. 2. Typical oscillatory stochastic dynamics of spins in 
the presence of anisotropie dissipation with a Drude bath, (a) 
A trajectory (solid) with parameters 7 = 400, B — 50/3T = 
50a;d, plotted on the Bloch sphere, (b) The same trajec¬ 
tory plotted with the unphysical choice of T = 0 to illustrate 
the deterministic part of the dynamics: The system oscillates 
around the 0(2) manifold with decaying amplitude. The lines 
0 = 0 and 9 = 9* (dashed) are shown, and an arrowhead in¬ 
dicates the initial state. 


Drude dynamics: The Drude bath has density of 
2 

states J{uj) = , where the Drude frequency, Wd, 

defines the bath cutoff. The new energy scale allows the 
bath to be in its classical limit, T 3> Wd, even when the 
equilibrium thermal fluctuations of the system are small, 
T B. It is then sensible to discuss the dissipation- 
dominated realaxation to this equilibrium distribution. 

The Drude bath has a memory on timescales deter¬ 
mined by the cut-off, ■ The dynamical dependence 
on the history of the bath is captured by introducing a 
time-dependent field term 

~ ^ (^^i ^diss) 1 

(18) 

Bdiss — (Bdiss [z • (s^ X Bj)] 7^) , 

where = 2^T6{t — t') and as before and we 

take ojd ^ T. This is shown in detail in Appendix B. 
The main effect of the bath is to induce oscillations in 0 
about the Markovian trajectory. In this strongly dissipa¬ 
tive limit the oscillations have a frequency and 

describe fluctuations that are small in the 9 direction and 
decay away on the bath memory timescale 1/wd. When 
the oscillation decay is much faster than the decay of 9 — 
i.e Wd 3> B j —we recover the Markovian 0(2) dynam¬ 
ics described by eqs (16). A derivation of this is shown 
in appendix B. The characteristic qualities of a typical 
trajectory are depicted in figure 2, whilst a simulation of 
the dynamics of equation (18) showing the separation of 
characteristic timescales is shown in figures 3 and 4. 

To conclude this discussion we summarise the phys¬ 
ical limits we have studied. We have assumed that the 
equilibrium thermal fluctuations are small B ^ T and 
that the environmental coupling is strong. We have fur- 
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(a) relaxation of 6 up to t = shown by the mean upper 

bound 9ij, mid-point 9m , and lower bound 9i^ of the trajectory 
envelopes. The the vertical line indicates the range of plot (b). 



^0 200 400 600 800 1000 


Bt 

(b) relaxation of </> up to t = /B shown by the mean upper 

bound ipu, mid-point and lower bound of the trajectory. 
(lOx rescaling) 

FIG. 3. Different timescales of relaxation in the stochastic 
dynamics of (a) 9 and (b) cj> when coupled to a Drude bath 
(Note the different plot ranges). An ensemble of 1000 spins 
(s = 1/2) initially at 0 = 37 r/ 4 , (f> = 37 r /4 when t — 0 
evolve with a magnetic field in the 9* = 7 r/ 4 , </* = 0 di¬ 
rection. The coupling is 7 = 5 x 10^, with energy scales 
B — lOT — lOOcUd, satisfying i? ^ T ^ Wd B/ys.On the 
timescales of 9 dynamics the fast oscillations in the trajecto¬ 
ries (see fig 4c) means that the trajectories are best charac¬ 
terised by an envelope with upper bound 6u, lower bound 9l 
and mid-point 9m, these are simply read off the oscillating 
trajectory as shown in fig 4c. Since the initial conditions is 
an extrema of the fast oscillations the initial point lies on 9u 
and 4>u. The ensemble averages {9m) (solid) and {9u), {9l) 
(dashed) are shown with {9u) +ou = {9u) + \/ {^u) ~ {^u)^ , 
and {9l) — ol = {9l) — \/ {dt) — (Sl)^ (both dot-dashed) 
illustrating the ensemble width, (a) The slow theta coordi¬ 
nate relaxes towards the equilibrium value 9* — 7 r /4 on a 
timescale 'ys/B, approaching it at t « 10'^ /B. The vertical 
line indicates the range of plot (b). (b) The same statistics 
are presented for the (p dynamics: The (p dynamics relaxes 
to its equilibrium distribution much faster on a characteristic 
timescale r ~ l/cud and is fully relaxed by t « 500/5. 


ther assumed that temperature prevents quantum cor¬ 
relations persisting in the bath, T 3> Wd-though this 
condition may be relaxed. Our analysis shows an effec¬ 
tive reduction in phase space from, 0(3) to 0(2), re¬ 
sulting from the confinement of the (p coordinate to a 
typically small region around (p = 0. This occurs when 
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(a) 9 relaxation on a long timescale, up to t = 10^/5. 
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(b) p relaxation on a shorter timescale, plotted up to t = 10^/5 
(lOx rescaling). 
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(c) very short timescale oscillatory behaviour in 9 induced by 
bath memory (see eq B5). 


FIG. 4. Different dynamical timescales of typical trajectories 
in 6 and tp when coupled to a Drude bath. A plot of three 
sample trajectories from the ensemble studied in fig 3. Being 
drawn from the ensemble these spins (s = 1 / 2 ) evolved in the 
same conditions: initially prepared aX 9 = 37 r/ 4 , (p = ‘iit jX at 
t = 0 and evolve with a magnetic field in the 9 = -n(A, cp = 0 
direction. The coupling is 7 = 5 x 10^, with energy scales 
B — IQT — lOOcJd, satisfying 5 ;§> T ^ tJd ^ Bj'ys. (a) 
Trajectories in 9 relax on a a long time scale. The vertical 
dashed line indicates the range of plot (b). (b) (p relaxes 

on a shorter timescale, the conhnement of (p is evidenced by 
the typically small excursions from (p = 0. (c) Oscillatory 

behaviour induced by the bath occurring on shorter timescales 
T ~ y/r^fuid is plotted for 9, similar behaviour occurs for 
(p. This behaviour is expanded on in appendix B, where the 
oscillations appear in eq B5. Each oscillatory trajectory can 
be characterised by the upper and lower edges 9u and 9l of 
its envelope, and its midpoint 9m. The ensemble statistics of 
these quantities are studied in fig 3. 
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B ^ T ^ (jjd ^ Bf^s, for both Markovian and Drude 
dynamics. In this parameter regime we have also shown 
that when coupled to the more general Drude bath, the 
system is described by the Markovian 0(2) dynamics of 
equation (16) on timescales greater than the bath mem¬ 
ory, I/Wd. 

V. CONSEQUENCES OF ANISOTROPIC 
DISSIPATION 

The Langevin equation derived in the previous sec¬ 
tion exhibits markedly distinct behaviours in different 
regimes. This has implications both for attempts to fit 
experimental data to such models and, ultimately, for the 
usefulness of a system described by them for computa¬ 
tion. 


A. Different behaviours for the same system 

For a qubit to be useful for quantum computation it 
must be sufficiently manipulable. For the simplest sys¬ 
tem of a single qubit, this in effect requires the imple¬ 
mentation of at least two non-parallel magnetic fields. 
The field, B, is thus assumed to be a tunable parameter 
of the system, whereas the environmental coupling, 7 , is 
fixed at some finite value. 

For the anisotropic coupling to the environment stud¬ 
ied here, the strength of dissipation depends not only 
upon the fixed parameter 7 , but also upon the orienta¬ 
tion of the magnetic field, B. This is important when 
characterising such a device. In particular a system that 
is analysed under conditions when B and z are nearly 
aligned will appear weakly coupled, whereas for other ori¬ 
entations of B the dynamics may be entirely dominated 
by environment-induced dissipative dynamics. 

It is a general feature of qubit systems that inhomo¬ 
geneous environmental couplings will result in dynamics 
that are correspondingly inhomogeneous^’^^’^®’^^. We see 
here that certain states evolve with minimal dissipation, 
whilst others are dominated by dissipative or noisy dy¬ 
namics. In extreme cases this may amount to a reduction 
in the effective state space. 

B. A model for lossy qubit arrays 

A natural application of our analysis is to understand 
some puzzling features of the D-wave machine"^^. This 
machine consists of a tuneable array of coupled Josephson 
junctions whose dynamics may be controlled to perform 
a quantum annealing or adiabatic computation. Vari¬ 
ous models have been posited for this system, includ¬ 
ing Bloch-Redfield simulations, Landau-Lifshitz-Gilbert 
models and two different 0 ( 2 ) models (a thorough review 
can be found in Ref 4). Our analysis sheds light upon the 
relationship between these different models and the fact 


that apparently quite different models yield surprisingly 
similar results. In particular, the salient features of the 
two 0 ( 2 ) models proposed in the under®- and over^’®^- 
damped regimes can be found in different limits of the 
Landau-Lifshitz-Gilbert model when proper attention is 
paid to the effects of the bath. These models are appro¬ 
priate in a limit where decoherence renders entanglement 
effects negligible. Bloch-Redfield simulations — the only 
one of the above to include entanglement — are expected 
to yield similar results in this limit. 

In the underdamped limit, the qubit precesses about 
the adiabatic minimum. Projecting this motion onto the 
polar angle results in harmonic oscillations of the polar 
angle about its adiabatic minimum. This is essentially 
the 0(2) model of Ref 5, though strictly the effective ki¬ 
netic term depends upon the local effective field felt by 
the quantum bit and so varys through the computation^®. 
This model is expected to accurately reflect the dynam¬ 
ics up to the point where it deviates markedly from adi¬ 
abatic, accounting for its success in predicting the prob¬ 
ability of correctly performing adiabatic computation in 
some circumstances'^®. 

In this manuscript we have focussed upon the over¬ 
damped limit. A model of over-damped classical rotors 
undergoing a thermal exploration of the 0 ( 2 ) state space 
was introduced in Ref 7. This model reproduced addi¬ 
tional, apparently quantum, effects®®, though it is inter¬ 
esting that the underdamped non-thermal 0 ( 2 ) model 
was also able to reproduce these results. Further statis¬ 
tics presented in evidence of quantum effects^ were repro¬ 
duced only by the over-damped classical model®^—and 
not the underdamped model. The overdamped model of 
Refs 7 and 24 can be obtained by artificially confining the 
motion of overdamped classical 0(3) spins to the 0(2) 
phase space^^. Our analysis in section IV B shows that, 
remarkably, anisotropic dissipation can bias the dynam¬ 
ics towards just such a confinement: the 0(3) dynamics 
of dissipative spin dynamics (Eqs. (11) and (18)) are re¬ 
duced to the effective 0(2) dynamics of Eq. (16) as illus¬ 
trated in Figs. 3 and 4. An unanticipated feature of this 
relaxation to an 0 ( 2 ) manifold is the restriction of the 
system dynamics to half of the 0(2) submanifold. The 
whole of this submanifold is a fixed point of the initial 
rapid decay, but half forms a stable and half an unstable 
fixed point. 

The representation of over-damped classical rotors 
in Refs 7 and 24 is at first glance rather different 
from the Langevin equation of Eq.(16). However, 
Metropolis-Hastings dynamics describes a dissipative re¬ 
laxation to thermal equilibrium and have been used pre¬ 
viously to model dynamics described by a Langevin equa- 
tion®®^"^®A5_ model of Refs 7 and 24 do 

not include the biasing of trajectories to half of the 0 ( 2 ) 
manifold that naturally arises from microscopic consid¬ 
erations. However, the dynamics of the component of 
qubit projected onto the direction of the local field is 
rather similar (see Appendix C), perhaps accounting for 
the success of Refs 7 and 24 despite their models not 



being strictly derivable from microscopic considerations. 

These analyses raise an immediate question of 
whether the D-wave system is in an overdamped or un¬ 
derdamped limit This is subtle. As discussed in V A, 
the strength of damping depends upon the microscopic 
details of coupling to the bath, the orientation of the ef¬ 
fective field relative to the z-axis, and the instantaneous 
position of the qubit on the Bloch sphere. Because of 
the latter effects, the dissipation is largest at the start of 
D-wave computation, when the effective field and Bloch 
spins are in the xy-plane, and decreases to zero as the 
computation proceeds. Even when the coupling to the 
bath is strong, and the initial dynamics overdamped, 
there is a transition to underdamped dynamics at some 
point in the computation. Whether or not the over- or 
under-damped dynamics determines the success or fail¬ 
ure of a computation depends upon precisely when this 
cross-over occurs. 

So is the D-wave system initially over-or under¬ 
damped? This question does not appear to be re¬ 
solved by published experimental data. To our knowl¬ 
edge no direct measurements of Ti or T 2 times on the 
D-Wave are available. Measurements of high frequency 
flux noise via macroscopic resonant tunnelling^® indi¬ 
cate that above a cutoff whf = 0.5 GHz noise is Ohmic 
with 7 hf = d)gS'$(0)/8/iA:BTL^ « 0.5, where the val¬ 
ues for inductance L = 265.4pH and shunt resistance 
R = 2fcrL^/S'$(0) = 20fcD measured in Ref 46, and is 
the magnetic flux quantum. However, our analysis shows 
that the qubit is also sensitive to noise with frequency 
lower than the system frequency. At the lowest frequen¬ 
cies, the noise is of 1// form^^. There is a large win¬ 
dow between the high frequency"^® and the low frequency 
measurements^^ in which the noise has not been directly 
measured. The measured current noise characteristics 
do not, therefore, preclude the possibility of 7 > 1 and 
initial over-damped dynamics. 


VI. CONCLUSION 

Our main result has been to show how anisotropic 
dissipation can bias quantum trajectories towards par¬ 
ticular sub-manifolds of the system’s Hilbert space. We 
have found a Langevin description of the dynamics of 
qubits that allows for anisotropic coupling to the envi¬ 


ronment. This is a natural generalisation of the Landau- 
Lifschitz-Gilbert equations which describe the dissipative 
dynamics of spins with isotropic coupling to the environ¬ 
ment. The fluctuation-dissipation relation has the impor¬ 
tant consequence that the anisotropic noise generated by 
this coupling inevitably leads to anisotropic dissipation. 

This model applies explicitly to qubits experiencing 
dissipation due to fluctuations in the level separation (en¬ 
vironmental coupling to the Sz operator). When the cou¬ 
pling to the bath is strong the anisotropic dissipation 
drives rapid relaxation to a reduced 0 ( 2 ) manifold of 
constrained dynamics. This emergence of this effective 
dynamics from the underlying microscopies reproduces 
some salient features of the dynamics of the models of 
Refs 7 and 24. These models were capable of reproduc¬ 
ing several observed behaviours of the D-Wave machine 
previously believed to evidence quantum dynamics. This 
highlights the necessity of understanding the dynamics in 
dissipative and strong coupling cases when interpreting 
the dynamics of an experimental system. 

Entanglement, which we neglect here, is crucial for 
full quantum dynamics, and necessary to get the ex¬ 
ponential scaling between the quantum state space and 
number of qubits. It has been argued recently in Refs 25 
and 26 to act as a resource for adiabatic computation. 
These works modelled quantum adiabatic computation 
with artificial constraints on the entanglement analogous 
to the artificial constraint of a local subsystem to an 0 ( 2 ) 
manifold. However one may anticipate that the effects of 
dissipation could naturally bias the trajectories to these 
restricted manifolds. Understanding this will be key to 
determining how best to use limited or dissipating entan¬ 
glement resources in computation. 

Understanding of the effects of state-dependent noise 
and anisotropic coupling to the environment is crucial 
for the proper control of quantum devices. As we have 
shown in the case of the D-Wave machine, these effects 
can bias the system dynamics in unexpected ways. Used 
constructively, this may be harnessed to useful ends. If 
ignored, the dynamics may completely different from that 
intended. 
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Appendix A: Strongly dissipative dynamics 


In this appendix we derive the effective 0(2) dynam¬ 
ics given in equations (16) and (17), in which the Marko¬ 
vian model features a confinement of </>. In the dissipative 
limit Tp ^ Td the spin dynamics are given by eq (14). To 
leading order in the long timescale, the (j) dynamics take 
the form 


(/) = —s")B sin 9* sin 9 sin (j) + rj = —^ sin (f) + rj. 

(Al) 

where as before = ysi? sin 0 sin 0*. This corresponds 
to a Fokker-Planck equation 


dp d /sin (j) m d \ 


(A2) 


for an ensemble distribution p. The equilibrium solution 
to eq (A2) is given by 


exp {A cos (j)) 
2 ttIo{A) 


(A3) 


where A = (yTr^)”^ and /„(•) are the modified Bessel 
functions of the first kind. Including the additional sub¬ 
leading terms (those in eq (14) that are missing from 
eq (Al)) leads to a Fokker-Planck equation with no closed 
form solution. However, the salient features are captured 
by making an appropriate shift to 4> and rescaling of A. 
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This leads to an equilibrium solution 


P = 


tan(^* = 


exp {Acos{(j) — (j)*)) 
27rIoiA) 
sin( 6 » - e*) 

57 sin^ 0 sin 0 * ’ 


where 9 is subject to a drift term originating from the 
mean value of the fast (p dynamics and a stochastic term 
which originates from the fluctuations in the (p dynamics 
away from their mean value. This approximation be¬ 
comes accurate on timescales t ^ Evaluating this 
exactly we obtain 


where A = (yTr^ cos which deviates from 

the exact solution only far from the distribution 
peak. In terms of 6, A takes the form A = 

sin^ 0 sin^ 0* + sin^ {9 - 9*) / (yTsind). 
Assuming a separation of timescales, the remaining 
9 dynamics can be found by averaging eq (14) over the 
equilibrium distribution of p given in eq (A4). This leaves 
a single equation giving the 9 dynamics, 

9 = —B sin 9* (sin p) -I- 77 ', (A5) 


(sin p) 


J dpp{p) sin 0 = sin( 5 i* 


hjA) 

loiA) 


(A 6 ) 


for the mean, whereas the noise term is defined by its 
mean = 0 and covariance 


WitWit')) 


{Bsin9*)'^ {{smp{t) sinp{t')) — {sinp{t)){sinp{t'))) 


{B sin 6 **)^ exp(—|t — t'\/T^) ((sin^ p) — {sinp)"^) 
2 t^{B sin 9*yS {t — t') ((sin^ p) — (sin p)"^') 


2T(j,{B sin 0 *)^(j {t 


t') 




hjA) 

AIoiA) 



+ cos^ p* 


hjA) - 
AIoiA)_ 


(A7) 


We consider the limiting cases in which the above 
form simplifies: For T B, we have A ^ 1 and 
Ii{A)/Io{A) = 1 — 1/(2A),-|-C)(A“^) which gives, to lead¬ 
ing order 


limit of T Wd, one obtains {'n{t)'q{t')) = T^{t — t') = 
Tywd exp (—Wd|t — t'\). This particular noise r]{t) can be 
written in terms of a (5-correlated stochastic dummy vari¬ 
able rj' as 


. M sin(6» - do) 

(sin0) = -d-; 

sj sin^ dsind* 

{v {t)ri'{t')) = 2'-fT{B t^ sin 9*)'^6{t - t'). 


(A 8 ) 


ri{t) = —Wd f dt'e 

J —00 


(Bl) 


For T ^ B, we have A < 1 and /i(A)//o(A) = A/2 -I- which satisfies {v'{t)ri'{t')) = 2'^T5{t - t'). Rewriting 
O(A^), which gives, to leading order in Td/rp, eq (8) the dynamics are obtained as 


(sin p) 


B sin(d — 6 * 0 ) 
2 T 7 sind 


(? 7 '(t)? 7 '(<')) = t^{B sin 9*)“^S{t - t') 


(A9) 


Substituting the values of eq (A 8 ) and (A9) into eq (A5) 
gives the forms eq (16) and (17) respectively in the main 
body of the paper. 


Si — Si X (Bi T Bdiss(t)) 

Bdiss(t) = -Wd f dt'e““‘^(*“‘^ (7Sj,2(t') - d') z- 
J —00 

(B2) 

Writing the second equation in its differential form we 
find 


Appendix B: Drude dynamics 


Bdiss — aid (Bdiss T d 

= -Wd (Bdiss - 7Z [z • (si X Bj)] - r/'z) 


(B3) 


In this appendix we derive the dynamical equation of 
a single spin coupled to a Drude bath and show that in 
an appropriate limit the long timescale d dynamics are 
given by the Markovian equation (16) with p remain¬ 
ing typically close to « 0. Using the Drude den¬ 
sity of states and evaluating eqs (9) and (10), in the 


where the second part follows by substituting Si^z and 
noting that z • (s^ x Bdiss) = 0. These dynamics equate 
to equation (18) after trivial relabelling. 

Separating eq (B3) into the fast {p, B^l^) and slow (d) 
degrees of freedom and solving the linearised equations 
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of motion for (j) and -B^iss 


6 = — B sin 9* (sin (p) +rj', 

/ • j,\ “ ^o) 

(sin0) =-^- 

S 7 sin^ 9 sin 9* 

{rj'{t)r]'{t')) ={B sin 6l*)^Cov(f, t') 


(B4) 


as before in equation (A8). However now there are 
long lasting oscillations in the noise, characterised by 
Cov{t,t') = (sin(/)(t) sin(()(t')) —(sin(/)(t))(sin(/)(t')). Eval¬ 
uating this one finds that 


=-fTT^{Bsin9*)‘^e 'x 

cos {uj\t — t'l) -I- — sin {uj\t — t'\) 
L 2w 


(B5) 


where a; = (l — characterises the oscillation 


frequency of the correlations. These fast oscillations can 
be seen in fig 4c. 

As expected on timescales much longer that l/wd the 
covariance has the same value as with the previous cases 


{ri'{t)ri'it')) k,2^T{Bt^ s\n9*Y5{t - t') 


2T 


7(ssind) 




(B6) 


thus on timescales much greater than the inverse bath 
frequency we recover the Markovian case (eqn (A8)) and 
the 9 dynamics of eqn (16) . 


Appendix C: Ensemble dynamics 

In this appendix we study the dynamics of an en¬ 
semble of non-interacting to show the appearance of 
the Markovian anisotropic dynamics at ensemble level. 
The ensemble dynamics are captured by a Fokker-Planck 


equation which is found to have similarly anisotropic dy¬ 
namics, which can in turn be related to Metropolis Hast¬ 
ings dynamics. 

We show that in the simplest case of a non-interacting 
spins the ensemble dynamics do not agree with the Bloch 
equations, this is clear as the anisotropic nature of the dy¬ 
namics persists even after the ensemble averaging. This 
result is perhaps surprising so we provide an explicit 
derivation from the Fokker-Planck equation. 

Within the first moment approximation the ensemble 
is described by the probability distribution 

P(s. t) = Z = y sinh (^s). (Cl) 

The first moment approximation in not appropriate when 
considering systems of interacting spins as correlations 
between the trajectories of different spins cannot be ne¬ 
glected. Thus for the purpose of deriving ensemble dy¬ 
namics, we consider a system of many spins, without in¬ 
teractions between them, each acting under the influence 
of an external field B and coupled anisotropically, along 
the z-direction only, to an ohmic bath. 

For a specific realisation of the history of the bath, the 
dynamics are given by the equation (11b). When we sum 
over the histories of the bath, the dynamics are described 
by the evolution of the probability distribution (Cl). 

Following the approach of Refs 16 and 19—in which 
an analogous calculation is performed for an isotropic 
bath coupling, i.e. microscopic dynamics corresponding 
to equation (lb)—we find that the evolution of p{s,t)is 
described by the Fokker-Planck equation 

p = Vs • [s X B -I- 7 S X z (z • [s X (B — TVs)]) ]p. (C2) 

By substituting equation (Cl) into equation (C2) and 
integrating over s, the dynamics of the parameter ^ 
are found to be governed by the equation 


i -I- I X B -I- 7s^ 




I X z [z • (I X B)] -I- yP 




(l - zz^) (I - lo) = 0. 


(C3) 


The fixed point of the dynamics is given by = B/P, 
which corresponds to the Boltzmann distribution, cr — 
(s) is the mean polarisation vector with norm = 

s ^coth (^s) — and cr' = ^ its derivative. 

The physical origins of the first three terms are clear, 
as they correspond to the respective terms of equa¬ 
tion (11b). The last term corresponds to longitudinal 


relaxation with a rate Ti = yP and transverse relaxation 
with a rate r 2 = It should be noted, however, 

that the dynamics described by these terms differs from 
the usual isotropic case due to the projecting out of the 
component in the z direction, this reflects the underlying 
anisotropy of the coupling. 

It is possible to re-write equation (C3) to give the dy¬ 
namics of the ensemble polarisation cr, where ^ is defined 
implicitly by the form of cr(^). Doing so, one obtains 
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<T + cr X B + 7 ^ ^1 - (T X z [z • (cr X B)] + 7 r (l - zzF') ^cr - = 0. (C4) 


Bt! y s 


0 12 3 4 
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FIG. 5. Comparison of Metropolis Hastings and anisotropic 
LLG dynamics. The ensemble averaged value (6) of 
Metropolis Hastings (dashed) dynamics is compared with the 
anistropic LLG dynamics (solid) data from fig 3. The LLG 
dynamics have small oscillations (see fig 4 for farther detail), 
thus we plot the ensemble averaged mid-point of these small 
oscillations, {0m), as previously in fig 3. Both simulation 
were performed using an ensemble of 1000 spins (s = 1 / 2 ) 
initialised a.t 0 = SnlA, evolving in a magnetic field in the 
0* — it/ 4:, (j>* — 0 direction, and a temperature T = B/10. 
For the LLG dynamics the coupling is 7 = 5x 10^, with energy 
scales B — lOT = lOOcJd, satisfying B ^ T ^ tUd ^ B/^s. 


For a highly polarised ensemble, u « s, at low tem¬ 
peratures, the ensemble dynamics converge upon the 
microscopic dynamics of equation (11b) showing that 
the unique behaviours described in the body of the 
manuscript persist in the ensemble dynamics of the sys¬ 
tem. These ensemble dynamics are shown in fig 3 where 
they are obtained from sampling the stochastic sampling 
of the trajectories defined by eq 11, a sample of these tra¬ 
jectories is shown in fig 4. In the main text we discuss the 
similarity between these dynamics and the Metropolis- 
Hastings dynamics of refs. 7 and 24. This similarity is 
most evident in their confinement to dynamics on an 0(2) 
manifold, but is clear also in the similarity between their 
dynamics, both of which are overdamped and dissipative, 
as shown in fig 5. The existence of such a relationship 
between overdamped spin dynamics and the Metropolis 
Hastings algorithm has been previously established^®“^°. 






